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The use of a relational time in quantum mechanics is a framework in which one promotes to 
quantum operators all variables in a system, and later chooses one of the variables to operate like 
a "clock". Conditional probabilities are computed for variables of the system to take certain val- 
ues when the "clock" specifies a certain time. This framework is attractive in contexts where the 
assumption of usual quantum mechanics of the existence of an external, perfectly classical clock, 
appears unnatural, as in quantum cosmology. Until recently, there were problems with such con- 
' structions in ordinary quantum mechanics with additional difficulties in the context of constrained 

' theories like general relativity. A scheme we recently introduced to consistently discretize general 

. relativity removed such obstacles. Since the clock is now an object subject to quantum fluctuations, 

OA ' the resulting evolution in the time is not exactly unitary and pure states decohere into mixed states. 

Here we work out in detail the type of decoherence generated, and we find it to be of Lindblad 
, type. This is attractive since it implies that one can have loss of coherence without violating the 

■ conservation of energy. We apply the framework to a simple cosmological model to illustrate how 

a quantitative estimate of the efi'ect could be computed. For most quantum systems it appears to 
be too small to be observed, although certain macroscopic quantum systems could in the future 
provide a testing ground for experimental observation. 
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I. INTRODUCTION 
A. Relational time in quantum mechanics 



(N 

' There exists interest in applying the rules of quantum mechanics in situations where there does not exist an obvious 
way of considering an external perfectly classical clock. An example is provided by quantum cosmology, presumably 
the correct description of the universe close to the big bang. Sometimes the discussion of quantum cosmology gets 
^H, further complicated by the problems associated with the quantization of general relativity. But it should be emphasized 
I ■ that the quantization remains troublesome even if one considers simplified model cosmologies, for instance described 
^jpj' by Newton's theory (an example is the Barbour-Bertotti model 

• • ■ Given the lack of an external time, one could try to use a variable internal to the system under study as a clock. Such 
. ^ ' ^ variable could not play the role of "t" in a Schrodinger equation, since one expects it to be on a similar footing as the 
other variables in the problem and therefore should be subject to quantum fluctuations. There have been proposals in 

^ the past to build a quantum mechanics using an internal variable as a clock. We will call these proposals "relational 
[ time" (see the next subsection for clarifications on terminology). What one does is to promote all variables in the 
theory to quantum operators, then choose from among the variables one (or several) that will operate as a "clock" 
and then compute conditional probabilities for the other variables to take certain values when the "clock" variable 
takes a given value. If one now considers the system as possessing an external perfect classical clock and builds an 
ordinary Schrodinger description for the system, and if the variable chosen as internal clock behaves semiclassically in 
such a description, then the relational picture and the Schrodinger picture will give similar descriptions of the system. 
If the variable chosen as "clock" has large quantum fluctuations, then both descriptions will differ. It is clear that 
since the "internal clock" will always have some quantum fluctuations, that both descriptions can never be identical. 
It is therefore of interest to address by how much both descriptions differ. This is one of the main purposes of this 
paper. 

In ordinary quantum mechanics, there is an old argument due to Pauli that one cannot promote a dynamical 
variable to a quantum operator associated with time, since this operator would have to be canonically conjugate to 
the Hamiltonian, and since time is expected to be monotonous, this would imply that the Hamiltonian cannot be 
bounded from below. We will see that in our approach this difficulty is not present. 

There is an additional complication. If one considers a system as alone in the universe, without external observers, 
there appear constraints among the variables of the system. Such systems are called "generally covariant" or "totally 
constrained" . For instance, in the Barbour-Bertotti model, which consists of a few Newtonian particles alone in 
the universe, one has the constraint that the total momentum, energy and the total angular momentum must have 
fixed values. In such systems, although one can formally write them as dependent on a parameter called "t" , one can 
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always reparameterize them and there is no natural choice of a variable to identify with a physical time. Sometimes 
this leads to calling them "reparametrization invariant" systems. In a cosmology based on general relativity, one 
has the Hamiltonian and momentum constraints. If one formulates the theory canonically in order to quantize it, 
the constraints have to have vanishing Poisson brackets with any quantity that can be considered physically. This 
presents problems, since it implies that these quantities are constants of the motion, and cannot be used as "clocks" . 
Page and Wootters 0| attempt to bypass this by choosing to build the relational framework in terms of quantities 
that do not have vanishing Poisson brackets with the constraints (that is, they choose to work at the kinematical 
level.) However, quantum mechanically, the states that are annihilated by the constraints are distributional within 
the space of kinematical states and do not lead to a good probabilistic interpretation. The resulting propagators are 
proportional to the Delta function and therefore "they don't propagate" as discussed in detail by Kuchaf Q. 

We therefore see that the presence of constraints interfere with the idea of attempting the introduction of a relational 
time in quantum mechanics. One possibility to circumvent the problem is to get rid of the constraints. We have recently 
proposed a discretization of field theories that has the property that when applied to totally constrained systems, the 
resulting discrete theories are constraint free. In spite of this, they manage in certain circumstances, to approximate 
well the continuum theories from which they were generated. This allows the introduction of a relational notion of 
time in the discrete theories. In this paper we would like to discuss, within the context of this approach, what are 
the quantitative effects of evolving a system using the relational time. 

It is clear that the relational time can be introduced in any context that would produce a system that is free of 
constraints. The discretized approach we consider is an example of such a context, but there could exist others. 



B. A note on terminology 

The idea of using a given dynamical variable in a closed system as a clock is quite old. One can find traces of it in 
the work of Leibnitz, Mach, Bergmann, DeWitt, Page, Wootters, Halliwell, Ashtekar, Rovelli, Smolin, Mermin and 
many others. For a partial list see 0,y, I3- The word "relational" is also used in an attractive proposal of Rovelli 
[Bi] in which observables (that is, quantities that have vanishing Poisson brackets with the constraints in constrained 
systems) are constructed via a relational technique. The proposal differs from the one presented in this paper. In it, 
the observables are constructed as dependent on a clock variable that is not promoted to a quantum operator like we 
do here. 



II. CONSISTENT DISCRETIZATIONS: BRIEF SUMMARY 



A. General framework 

We summarize here the consistent discretization approach, for a lengthier discussion see Readers whose main 
interest is not the details of how to handle general relativity can skip this section. We start by considering a 
system with an action as in ordinary classical mechanics dependent on a finite number of configuration variables q"', 
S = ^ dTL{q°' ,q°'). Although this may appear oversimplificatory, recall that when one discretizes a field theory one 
is really dealing with a mechanical system. This Lagrangian should be considered as a very general Lagrangian, in 
the sense that it can accommodate first order formulations (in which case the momenta are considered as g's) and 
systems with explicit constraints in the Lagrangian (in which one treats the Lagrange multipliers as g's such that 
no q appears in the Lagrangian). The formalism can accommodate theories invariant under reparameterizations, like 
general relativity, for that reason we refer to the evolution parameter as r. 

We now discretize the evolution parameter (in the case of a field theory we also discretize space), and replace the 
derivatives via if^ = {q'^j^i — g^J/Ar (other discretizations of the derivatives are possible, but would require reworking 

the formulation). The action will then become a sum S = X^iLo ^{qn, q-n+i) where we have absorbed the interval of 
discretization into the definition of L{q,qn+i). The Lagrange equations are obtained by minimizing the action with 
respect to each of the , 

9L(g°,g°^i) dL{q-_,,q-) _ 

We next introduce a type 1 canonical transformation that corresponds to the evolution of the system from n to 
n + 1 with generating functional —L{qf^, Qn+i)^ 
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It is immediate to see, by evaluating the top equation ai n — n ^ 1 that the two equations are equivalent to the 
Lagrange equations. These two equations should be seen as an imphcit map between and qn+i^Pn+i- If the 

implicit map is invertible, then the resulting transformation is indeed canonical, by construction. 

If the map is not invertible, that is, if the determinant det g^^J^ — vanishes identically, one has constraints. We 

denote them generically as (/)'^((7j^,p°), yl = 1 . . . M. The variables at instant n+1 are not completely determined by 
the variables at instant n, and will depend on some free parameters , B ^ 1 . . . A/, 

qn+l = qn+l{qn,Pn,V^) (4) 
= Pn+l{qn,Pn,V^). (5) 

It should be noted that the constraints not only occur at level n but at all levels. We will treat them by considering 
the constraints at level n and then ensuring conservation upon evolution. This can lead to four different situations: 
1) imposing the constraints at all levels leads to the determination of some of the indeterminate parameters V^; 2) 
new (secondary) constraints appear; 3) the constraints are preserved automatically; 4) the system is incompatible. In 
the three compatible situations one generically ends up with evolution equations, potentially depending on some free 
parameters and potentially with extra constraints. It can be shown that the resulting evolution equations preserve 
the Poisson brackets and (weakly) the constraint surface. 

It should be noted that preserving constraints in the discrete theory is harder to accomplish than in the continuum 
theory. In the latter one only needs to show conservation infinitesimally, and the calculations neglect higher order 
terms in the (infinitesimal) time interval. In the discrete theory, since the evolution is finite, one needs to show 
conservation preserving all powers of the time interval. For this reason, in many circumstances constraints that are 
first class in the continuum become second class when discretized. In our approach, constraints are eliminated by 
solving them for the Lagrange multipliers. One therefore ends up with a discrete theory that is free of constraints 
that nevertheless approximates a continuum theory that has constraints. In the next subsection we illustrate this 
procedure in detail by considering an example in general relativity. 

To quantize the theory one has to deal with the constraints that may be left, since as we mentioned it may not be 
always possible to solve the constraints by choosing the Lagrange multipliers. If the remaining constraints are second 
class one imposes them strongly and one replaces Poisson brackets by Dirac brackets, if they are first class they are 
imposed as operatorial identities on the wavcfunctions. Finally, one needs to implement the canonical transformation 
as a unitary evolution operator. 



B. A cosmological example 

We will consider a Friedmann cosmological model, written in terms of Ashtekar's variables 0. The fundamental 
canonical pair is {E, A) where E is the only remnant of the triad after the minisuperspace reduction and A is its 
canonically conjugate variable. We will consider the presence of a cosmological constant and of a scalar field. We will 
assume the scalar field has a very large mass so we can neglect its kinetic term in the Hamiltonian constraint, for the 
sake of computational simplicity. The Lagrangian for the model is, 

L^EA + TT4)^NE'^{-A^ + {A + m^(f>^)\E\) (6) 

where A is the cosmological constant, m is the mass of the scalar field (f>, ir is its canonically conjugate momentum 
and N is the lapse with density weight minus one. The appearance of \E\ in the Lagrangian is due to the fact that 
the term cubic in E is supposed to represent the spatial volume and therefore should be positive definite. 
We consider the evolution parameter to be a discrete variable. Then the Lagrangian becomes 

L{n, n + l) = En{An+i - A„) + ^„((^„+i - 0„) - NnEl{-Al + (A + m^cbl)\En\) (7) 

If one carries out the construction of the previous subsection, one ends up with the following evolution equations, 
(see 111 for details), 

F„^+i = 4^,6-1 (8) 
An+^ = ^^^7^ (9) 

4>n+l = 4>n (10) 
Pt+1 = Prt - K - Pn^^n) m'cf>^Q~^ (11) 
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Q = A + m?(j)'^. The variable E and its canonical momentum have been eliminated, and we have solved the constraints 
of the model for the Lagrange multiplier N. This evolution preserves the Poisson brackets among the canonical pairs 
iA,P^) and (0,P'*). 

The classical solution of the finite difference equations can be seen to approximate well the continuum solution (see 
for details). The approximation gets better for later times, and one enters a "continuous regime" in which one 
can approximate the dynamics by a differential equation in terms of a continuous variable. Such equation can be 
integrated exactly and from there we can read the behavior of the variables asymptotically into the future. The result 
is, 

A - Wck\/Aa(n + fc)2/3, (12) 

E = 4la„ck«'(" + fc)'/', (13) 

where a and k are integration constants, which we chose to be dimcnsionless and therefore we made explicit the appro- 
priate dimensions through the introduction of the Planck length ^pianck- The connection A is therefore dimensionless 
since the cosmological constant A has dimensions of (length) 

The model can be quantized by implementing the canonical evolution equations as operator equations via a unitary 
evolution operator, 



< A', 0'|C/|A, sg(A)^^e-s(-4)A^^^(^/ _ ^) (14) 

Dirac @ had already noted in 1933 that the unitary operator that implements a canonical transformation is given 
by exp(— iG) where G is the generating function of the canonical transformation. In our case the generating function 
(after eliminating the Lagrange multiphers is indeed given by G{An, An+i) = A^{An+i — An)/Q. There is an overall 
difference with Dirac's result since he chooses a specific factor ordering that does not coincide with the one we chose. 
It is interesting that this construction is what led Dirac to the notion of path integral. 

To define a time we introduce the conditional probabilities, "probability that a given variable have a certain value 
when the variable chosen as time takes a given value" . For instance, taking A as our time variable, let us work 
out first the probability that the scalar field conjugate momentum be in the range AP"^ = [P^"^^ , Pj-j-)] and "time" is 
in the range AA — [^(i), A(2)] (the need to work with ranges is because we are dealing with continuous variables). 
Since we have no constraints, the wavefunctions ^'[A, 0, n] in the Schrodinger representation admit a probabilistic 
interpretation. Therefore the probability of simultaneous measurement is, 

Psim(AP'^,A^) - lim — y / ^'^[A,P'^,n]dP'l'dA. (15) 

We have summed over n since there is no information about the "level" of the discrete theory at which the measurement 
is performed, since n is just a parameter with no physical meaning. With the normalizations chosen if the integral in 
P"^ and A were in the range (—00,00), Ps\xa would be equal to one. 

To get the conditional probability Pcond(AP'^| A^), that is, the probability that having observed A in A^ we also 
observe P"^ in AP"^, we use the standard probabilistic identity 

P,i„,(AP^ AA) = P(AA)Peond(AP^|AA) (16) 

where P(AA) is obtained from expression (|15|l taking the integral on P"^ from (—00,00). We therefore get 

Pcond(AP'^IA^) = '^^^ . (17) 

Notice that all the integrals are well defined and the resulting quantity behaves as a probability in the sense that 
integrating from (—00,00) in P"^ one gets unity. 

We will return to this example when we attempt to estimate the effects of decoherence in our current universe. 



III. RELATIONAL TIME IN GENERALLY COVARIANT SYSTEMS 



Let us return now to the general discussion of the use of a relational time in generally covariant systems. We 
assume we have discretized the system of interest using the techniques describe in the previous section. Therefore the 
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resulting discrete theory is either free of constraints, or if, there are constraints left, we need to ensure we work with 
physical quantitites that have vanishing Poisson brackets with the remaining constraints. An example of this would 
be general relativity written in terms of Asthekar variables in a generic situation. There our technique eliminates 
the diffeomorphism and Hamiltonian constraint but leaves the Gauss law as a constraint. One therefore would have 
to build the arguments that follow involving quantities that are invariant under gauge transformations, for instance 
using Wilson loops. 

As we discussed in the previous section, we have a unitary operator to be related to the canonical transformation 
on the phase space qmPn <ln+i,Pn+i- We shall denote this operator as U{n,n') for a general displacement n n' . 
At this point it is worthwhile re-examining the Pauli argument against promoting time to a quantum operator. Since 
in the discrete approach the evolution of the system from n to n + 1 is done through a unitary transformation, 
corresponding to a discrete canonical transformation in the classical theory, the evolution is not associated to a 
Hamiltonian. One can define (locally in time) objects that behave close to a Hamiltonian at least when one considers 
configurations that approximate the continuuum well. But such objects do not exist globally |lC| . Therefore one 
generically does not have a Hamiltonian with which to build the argument laid out by Pauli, which is of global nature. 
Notice that in terms of the parameter n the evolution is globally defined, but n is not what is promoted to an operator 
representing time. 

Let us introduce now the set of self-adjoint operators On{qmPn) with eigenvalues o on the isomorphic set of Hilbert 
spaces Hn — L?[qn)- We can construct now the set of projectors Po{n) = J^^'^j \o,j,n >< o,j,n\, where we have 

assume a continuum spectrum, and denote the eigenvalues of all other operators that form a complete set with O, as 
J- 

We will construct the conditional probability to measure the value o (since we consider a continuous spectrum, one 
should strictly speak of measuring a value within an interval Ao surrounding the value o) for the partial observable 
0{qn,Pn)- We also consider a "clock variable" T{q,p) and the associated self-adjoint operator Tn{qn,Pn) with eigen- 
values t. We define the projector Pt{n) = ^^^^f^ \t,k,n >< t,k,n\ as we did for 0„. In this case k denotes the 
eigenvalues of all the operators that form a complete set with T. These are well defined projectors for each n value, 
that is Pa{nf = Pa{n),Pain)Pa<in) = S{a - a')Pa{n),Y.a PaH = 1 V n. 

We are now ready to introduce the relational time. For an application of these ideas in the simple example of the 
parameterized non-relativistic particle see [T]| . The relational interpretation is defined as follows. Let us assume that 
the system is initially in the state described by the density matrix p. The conditional probability to obtain the value 
o E Ao for the quantity O given the value t E At for T is. 



E„ Tr(Po(n)P,(n)pP ,(n)) 



V{o e Ao\t e At)^ = ^^'^Zlrr^/ZV - (18) 



It should be emphasized the important role that the presence of the parameter n plays in these formulas. Without 
such an ordering parameter, one could not define the conditional probabilities. The parameter n introduces a notion 
of simultaneity in the construction (at a given spatial point). In fact, several previous attempts to introduce relational 
times were problematic due to the lack of such a parameter. For instance, this led Unruh jlJl to attempt to introduce 
in an ad-hoc manner a "mysterious time" in the continuuum theory to play such a role. 

If we further assume that O and T commute, we can construct the projector 



PoA^)-J2 I 

I 

and rewrite equation (|18|l as. 



\o,t,n,l X o,t,l,n\ (19) 

At 



7^(oeAo|^eAt)^^ ^^"f/f^(")^) (20) 

^ E„ Tr( j doPo,t(n)p) 

and from now on for simplicity we will assume that O and T commute. 

Due to the unitary evolution and the cyclic property of the trace, if we define the operators Ho = Po(ri),Ht = 
X^n Pt{ii') a-nd Ho t = Pojtin), we can rewrite the conditional probability as. 

The reduction postulate is given by, 

P^$^^^o,t(ri)pPo.t(ri). (22) 
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This reduction process allows us to calculate the conditional probability that defines the correlation functions 
(propagators), 

'^^'^ ''^'^'^^ = E„,.Tr(P.(n')P.,(T),P„,(n)) ' ^''^ 

The difference with usual propagators is that here the times t,t' are the outcome of a quantum measurement. 
Therefore this will lead to a meaningful definition of probability if we interpret the above expression as having 
measured o,t and t' and asking what is the probability that one will measure o' given those measurements. Notice 
that at the moment there is no well defined concept of "time ordering" , since the latter is only expected to arise in 
semiclassical regimes. There is however a well defined notion of simultaneity (at a given spatial point). In this sense 
the measurements of o and t should be simultaneous and so should be the measurement of t' and o'. 

In non-relativistic quantum mechanics one can compute the probability that measurements of position will find the 
particle at a position x at a series of times ti,t2, ....t„. Here one can compute such probabilities through a natural 
extension of (|23|l . This introduces a reduction postulate in the framework and leads to the history approach for 
computing a probability for a succession of events. This works if one chooses as clock a robust variable, i.e., such that 
the information about the variable is not destroyed in the measurements. 

From the previous expression, we have that if one prepares a quantum state with eigenvalue o at time t, po.t the 
probability of it to evolve into a state with eigenvalue o' at time t' can be written as. 



Since up to now we have made almost no assumptions about the nature of the "time" t chosen, it could happen that 
the same value of t occurs many times upon "evolution" in the parameter n. This eliminates the predictive power of 
the theory, at least locally in time in the following sense: one could make a definite prediction only upon completing 
the entire evolution of the system and determining if the variable t takes a given value more than once. Only then 
one could make sense of the probabilities and predict the probability of a given observable taking a given value at 
"time t" . We will discuss this in detail in a forthcoming paper. 

As we shall see in what follows, there is a particular regime where we can concentrate on a particular range of steps. 
This regime will correspond to considering as a clock a variable that operates semiclassically. When this happens, then 
the quantum mechanics we constructed will reproduce the results of ordinary quantum mechanics. The agreement 
will be better the more classically the clock variable behaves. 

Another element to be emphasized is that we have used the parameter n to define simultaneity. One has to be 
careful in the case of general relativity that one has discretized the theory in such a way that the parameter n is 
associated with Cauchy surfaces of the continuum theory. Otherwise it would not be correct to use the parameter to 
define a notion of simultaneity. Even this requirement is not enough. In the consistent discretization scheme the lapse 
is determined dynamically. One cannot rule out situations in which the lapse becomes negative. In such situations 
"time runs backwards" and one covers the same region of the manifold more than once. These situations in the 
classical theory should also be avoided in order to obtain a sensible quantum mechanics. It should also be emphasized 
that the whole issue of the covariance of the formalism is still to be worked out in detail. In particular, since we are 
here largely concentrating in models without spatial degrees of freedom, it is difficult to even pose the question of how 
they transform under coordinate transformations. We are taking into account possible time reparameterizations, but 
not coordinate changes that mix space and time. We are currently studying model systems in which some limited set 
of symmetries can be implemented and study how the relational approach will mix with covariance, but we will not 
discuss this in this paper. 



IV. SEMICLASSICAL TIME 



We now wish to show that the new quantum mechanics we have created reduces to ordinary quantum mechanics if 
the time variable chosen behaves classically. If one takes into account semiclassical corrections for the time variable, 
then one will end up with corrections to ordinary quantum mechanics. At this point it is worthwhile mentioning the 
work of Egusquiza, Garay and Raya [13| . They have studied modifications to quantum mechanics through the use of 
imperfect clocks. They treat the clocks classically but admit they may have fluctuations in their behavior, perhaps of 
quantum mechanical origin, and they model them through a Markovian process with a given probability distribution. 
There are many points of contact between their calculations and the ones we present here. The main difference relies 
on the fact that we are treating the clock as quantum mechanical and in particular we allow its probability to evolve 
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as a function of time. We are not including any thermal or other effect, but these could definitely be incorporated if 
needed through the formalism of Egusquiza et al. 

Let us assume now the existence of a semiclassical regime for the variable chosen as time for a given initial state 
of the complete system p. That means that there exists a region R in the spectrum of the operators T„ such that for 
a value t € R there is an interval of values Atn of the step parameter n such that Vn{t) ~ OVn ^ Atn. Also given a 
value of the step parameter no, there exists an interval around it A„pt such that Vnoit) ~ V f ^ ^no^ = 0. In these 
expressions Vn{t) is the probability density of having the value t at a given n. We speak of probability density since 
we are assuming that the operator T has a continuous spectrum. If the spectrum were discrete then Vn{t) would 
be the probability of having the value t at the level n. This semiclassical limit therefore implies a strong correlation 
between the n parameter and the eigenvalues of the observable T in region R. Notice that in this case each reduction 
process is almost a projection given by, 

Po.M)pPo.M)- (25) 

We need to assume also that the clock and the system in study are weakly interacting (otherwise one would not 
have a reasonable clock). Concretely, we assume that the state of the system is of the form of a product p ~ pi 0/92 
throughout the evolution determined by a unitary operator U « 0) f/2 also of product form. We take "1" to 
describe the clock and "2" the physical system under study and the weak interaction also implies that Ui and U2 
commute. 

Up to now we have considered the quantum states as described by a density matrix at a given level of the dis- 
cretization parameter n. Since the latter is not an observable, we would like to transition to a description where we 
have density matrices function of the observable time instead of n. In order to do this, let us recall the expression 
for the usual probability in the Schrodinger representation of measuring the value o at time t in ordinary quantum 
mechanics, 

^^"'^)^= Tr(p(^)) ' ('^^ 

where the projector is evaluated at t = since in the Schrodinger representation the operators do not evolve. We 
would like to obtain an expression similar to this one in the relational time picture. We are denoting p[t) as the 
density matrix in the picture in which we have a time variable, we will drop the tilde to denote the density matrix 
that arises under the evolution of the step parameter n and we will drop the n dependence to refer to the initial 
density matrix p(n = 0). 

We start by considering the conditional probability defined in (|20|l . 

}J„Tr(Pt(n)p) 

and one could have omitted the last Ptin) using the cyclicity of the trace since we are assuming that P^ and Pq 
commute. Given that we are assuming the presence of continuous spectra, the projectors in the above expression 
should be understood as integrated over the interval, i.e. Po{n) = J^^Po'{n)do' and similar for Pt. We now introduce 
the hypothesis that the clock and the rest of the system interact weakly and write explicitly the evolution of the 
projectors in the step parameter n to get. 



V{oG Ao\t e At)p = 



En Tr (ul{n)Po{0)U2{n)Ul{n)Pt{0)Ui{n)pi ® p2 



E„Tr(Pi(n)pi)Tr(p2) 
En Tr (ul{n)Po{0)U2{n)p2) Tr (ul{n)Pt{0)Ui{n)p, 
E„Tr(P,(n)pi)Tr(p2) 



(28) 



From this expression, using the cyclic property of the trace, we can identify the expressions of the density matrix 
evolved in relational time. We start by defining the probability that the measurement t corresponds to the value n, 

TT(Pt{0)Ui{n)piUl{n)) 
E„Tr(Pt(n)pi) 

and notice that En'^ri{t) = 1- 
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We now define the evolution of the density matrix, 

~P2{t) = U2{n)p2Ulin)Vn{t), (30) 

n 

and noting that 

Tr (p2(<)) = ^n.(t)Tr (p^) - Tr(p2) (31) 

n 

one can equate the conditional probability l|28|) with the usual expression for a probability in quantum mechanics 
(|26|l . It should be noted that all the sums in n, due to the assumption that the time variable is semiclassical are 
only nontrivial in the interval A^n since outside of it, probabilities vanish. Something else to notice is that when we 
introduced the projectors, there was an integral over an interval. Therefore in the above expression for the evolution 
of the density matrix, this has to be taken into account. Since the interval At is arbitrary, one can consider the limit 
in which its width tends to zero, apply the mean value theorem in the integrals, and the interval in the numerator 
and denominator cancel out, yielding an expression for p2(0 that is independent of the interval, and involves the 
non-integrated projector -Pt(O). 

We have therefore ended with the standard probability expression with an "effective" density matrix in the 
Schrodinger picture given by P2(t). In its definition, it is evident that unitarity is lost, since one ends up with a 
statistical mixture of states associated with different n's. We also notice that probabilities are conserved, as can be 
seen by taking (|26|l and integrating over x. We recall that p2 is not the normalized density matrix; the latter can be 
easily recovered dividing by the trace. 

We will assume that Vn(t) = f(t — tmax{n)) with / a function that decays quite rapidly for values of t distant of 
the maximum tmax which depends on n. 

To manipulate expression (|30|l more clearly, we will assume we are considering a finite region of evolution and we 
arc in the limit in which the number of steps in that region is very large. We denote the interval in the step variable n 
as going from zero to N with N a very large number. We define a new variable v — en with dimensions of time such 
that Ne — V with V a chosen finite value. We can then approximate expression H30fl by a continuous expression, 

P2{t)= I dvf{t-traax{v))p2{v). (32) 

In this expression imax(i') = traa.^{n = v / e) and 

P2{.v) = [/2(n - v/e)p2Ul{n = v/e). (33) 

In all the above expressions, when we equate n — v /e it should be understood as n = Int(i;/e), which coincide in the 
continuum limit. (Notice that strictly speaking we should write p2{v /e) to keep the same functional form as for P2{n), 
but we will drop the e to simplify the notation.) 

Let us assume now that tmaxiv) = v + eTi{v) + e^r2(v), with ri(D) ~ 1 and r2(t;) 1, that is, the value at which 
Vn {t) has its maximum depends linearly on n plus a small correction. One can always obtain the classical dependence 
of t on n solving the equations of the theory. The linear approximation can be justified locally as a Taylor expansion, 
or by redefining the variable chosen as time such that the relation with n is linear. 

To simplify calculations we replace the function / with a function / which agrees with it up to terms of order e^, 

fit - W(«)) = fit~v- eTiiv) - e^r2{v)) = f{v - t + eV^it) + e'r2it) - e^aTiit)), (34) 

where a{v) = dTi{v)/dv and we have neglected terms of order e"^. 
We now assume that / is well approximated by a Dirac delta, 

J[v - t + eri(i) + €'T2{t) - e'aTiit)) = co{t)5{v - t + eTi{t) + (^T2{t) - e^aTi{t)) 

+a(<)e(5'(w - t + eVi{t)) + h{t)e^&"{v - t) + 0{(!^) (35) 

with 6 > 0. Since Vvit) = f{v-t + eTi{t) + e'^T2it) ~ e'^aTi{t)) and since we want Vv{t) = 1 then co(t) = 1. 
We now substitute the explicit form of / H35|l in 132|l . which allows to compute the integral in v explicitly, 

o o2 

hit) = P2{t - eV^it) - e^T2{t) + e^a{t)Ti{t)) - a{t)e—p2{t ~ €T^{t)) + b{t)e^—p2{t) + 0{e^). (36) 
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We proceed to expand the arguments of p2 in this expression in powers of e. We also use the fact that in the 
continuum Umit we can now associate a Hamiltonian with the evolution operator U2(v) = eKp{iH2v), and recalling 
that the evolution of the density matrix is P2{v) = U2{v)p2Ul{v), we get, 

hit) = P2{t) + (leir^it) + a{t)))[H2, P2{t)] - le^ («(i)ri(t) - Ta^) [H2, P2it)] 

-e\^T,itf + a{t)T^it)+bmH2,[H2,P2im. (37) 

Notice that generically, one can write the discrete generator as U2{n) = exp{iH2^^n). However, there may be points 
where the canonical transformation is singular, and example is the Big Bang in the cosmological model |lOl| . the 
series defining the logarithnm of U2 fails to converge and therefore H2 does not exist. Notice also that in taking the 
continuum limit TJ^'^/e = H2. 

We will find useful to have the inverse relation between the density matrices, 

P2it) = Mt)-{ieiri{t) + am[H2,P2it)]-eH^iit)+ait)f[H2,[H2,P2im (38) 
(a(i)ri(i) - r2{t)) [H2,p2{t)] + {^^litf + ci{t)T,{t) + b{t)^ [H2, [H2, PslO]]- 

We now take the time derivative of equation (|37|) and we note that the derivatives will act on the coefficients 
a,b,a,T{t) or on P2{t). We replace P2{t) with the commutator with the Hamiltonian. The final result, using (|38|1 to 
rewrite everything in terms of p2 is 

9p2{t) ( d ,^ , ^ , 25(a(i)ri(t) -r2(t))\ ^ , Adfaitf ,,\\,,, . , 

^^^i[-l + e-{T,{t)+a{t))-e^ dt ) [^2, P2(t)] + ( 2 "^^^VJ t^^, [^^2, P2(i)]]- 

(39) 

The general form of the resulting evolution equation for the system under study is therefore, 

= -l[{l+m)H2rP2] - <y{t)[H2, [i?2,P2]], (40) 

with the coefficients (3{t) and a{t) that are functions that are small corrections (in terms of the expansion in e), and we 
have neglected terms involving triple and higher commutators with the Hamiltonian since they correspond to higher 
powers of e 

The presence of (3{t) just redefines the Hamiltonian of the theory. The coefficient a{t) is the one that gives rise 
to new effects. In particular it will imply that pure states evolve into mixed states. In fact, the equation we have 
obtained is a particular case of equations that are considered in the context of decoherence in quantum mechanics, 
called Lindblad type equations 

^^=.^[H,p\-V{p), (41) 

with X'(p) satisfying the properties, 

I?(p) =^[D„,[D„,p]], D^^Dl [i?„,i/] = 0, (42) 

n 

SO it defines a completely positive map on p that is consistent with the monotonous increase of Von Neumann entropy 
S = Tr(plogp) and conservation of energy. This type of equation was introduced by Ghirardi, Rimini and Weber 
(GRw) m with the aim of providing an objective solution to the measurement problem in standard quantum 
mechanics. (Similar equations can be used to describe the decoherence due to interaction with an environment, see 
Ra|.) GRW considered a single Z? as a localizing operator in coordinate space. As discussed by Adler and Horwitz 
|17|. and also Milburn, Percival and Hughston IT^. setting D to be proportional to H is most natural since it 
leads to an objective state vector reduction in the energy pointer basis. This loss of coherence may be a way to avoid 
macroscopic superpositions, like the "Schrodinger cat" 17, 18]. 

We can recognize that in our case there is only one Dn that is given by the Hamiltonian. Having a Lindblad form 
is desirable since it implies that the Hamiltonian will be conserved automatically by the evolution considered. Other 
proposals for decoherence from quantum gravity, like Hawking's $-matrix may have problems with the conservation 
of energy poj. 

In order to study the influence of (3{t) and cr(t), we start by assuming that they are constant as functions of time 
f3{t) = f3 a{t) = a. In general these functions will be constant will some small fluctuations, in which case we can 
interpret the bars as average values. In such case the equation can be solved exactly, 

P2nmit) = P2„„(0)e-(i+^)--*e(-*(--)^)*, (43) 
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where we have written the density matrix element in the basis of energy eigenstates, and Unm is the Bohr transition 
frequency between the state n and m in the energy basis. It is important to notice that a has to be positive for the 
evolution to be physically acceptable. Otherwise the trace of the square of the density matrix will be larger than one. 
Let us estimate the sign of 

ait) = 

The coefficient a(t) represents the asymmetry of the probability of the time t as a function of the evolution parameter 
n. There is no reason for this probability distribution to have a definitive asymmetry, so its mean value through the 
evolution will vanish. The coefficient hit) represents the spread of the probability distribution. Suppose one starts 
the clock in a quantum state that is peaked around a certain value tg. As the parameter n evolves, the quantum state 
will disperse and the spread of the probability of f as a function of n will disperse too. Therefore 6(t) will increase 
and its derivative will be positive. The magnitude of the effect will depend on the details of the clock. For instance, 
if one used a free particle as clock the spread of the wavepacket will be linear and therefore a(t) will be a constant. 
The important point is that it is positive, and therefore signals an arrow of time. 

In the real world, the subsystem chosen as a clock will be subject to several other kinds of fluctuations, for instance 
thermal fluctuations. The influence of these kinds of errors leads 111 to an equation similar to H40|) . but where the a{t) 
is proportional to o? . We therefore see that a quantum evolution will exhibit decoherence from fundamental effects, 
like the ones considered in this paper, and also environmental effects. The latter could in some cases be minimized 
by choosing the experimental setup in appropriate ways, the fundamental effects will, however, be an ultimate limit 
to the achievable coherence in a quantum evolution. 

It should be noted that the value of a depends the choice of clock, and on the quantum state chosen for the clock 
initially, and on the interaction of the clock with the environment. 



d_ fa{ty 
dt 



bit) 



(44) 



V. ESTIMATING THE MAGNITUDE OF THE DECOHERENCE 



We will now proceed to use the cosmology we discussed earlier to give an estimate for the level of decoherence. This 
can only be viewed as a first calculation, a more realistic model would be desirable. In particular, we have chosen a 
model with a scalar field in the infinite mass limit, in which the model is a DeSitter universe. This is not a very good 
model for approximating the universe since it is homogeneous in time (and in particular the "Big Bang" we refer to 
is just a coordinate singularity). One should think of it only as asymptotically approximating our universe into the 
future. It should only be taken as a guide of the way to compute a. The calculations could be repeated in a more 
realistic model, but would be more involved. 

The task at hand is to provide a justification for the lattice spacing of the discrete model. It is evident that 
such spacing could not be determined if we just considered a non-gravitational system as our only object of study. 
Considering the influence of gravity will allow us to introduce a fundamental length scale in the problem, the Planck 
scale. 

Recalling the deflnition of b in l|35|) we can relate it to the quadratic deviation of the time variable from its linear 
behavior t — en, due to the spread as a function of time of the wavepacket representing the clock through 



therefore we can estimate a to be, 



.^i,^i^,JAl!l„ (46) 



(we will neglect factors of order one in working out the estimates) . 

We now need to choose a time variable in the cosmology. As before, we will choose it related to the connection A, 
but we also wish to have a linear relation between time and n, as we assumed in the calculation of the decoherence. 
Therefore we rescale 



t= A^l\ (47) 



A3/4/'l/2 

*Planck 



which taking into account the asymptotic behavior of the cosmological variables we described in section (jllBII yields, 

t = a3/2^Pianck(n + A:) (48) 
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from which we can therefore read off e = a^^'^ipia.nck and therefore Ae = |-v/aAa^pianck- We are choosing the 
integration constant k to be very small compared to n, since we are assuming we are in the asymptotic future region. 
It should be noted that the singularity happens for n + fc '--^ 0. 
So we can therefore estimate u, 



(Aa) 



a = ^'^"'^'^ ^ '-n. (49) 

e 

We can now proceed to find a bound on a(Aa)^ based on the uncertainty principle, AEAA > ^lian^k' which yields, 

a(Aa)2 > (50) 

SO we finally get for the decoherence parameter, 

^Planck ^Planck / r -i \ 

a > — -=— = — (51) 
eVAn VAt 

At this point it is worthwhile pondering what is meant in this equation by "t" . We have identified this variable as 
related linearly to the parameter n via equation (|47|l . This choice involves selecting a prefactor with dimensions of 
time. This prefactor is arbitrary. Therefore it would appear that the result for a is dependent on this factor. Indeed 
it is. However, it should be noted that a similar factor arises in identifying the continuum limit Hamiltonian H2 from 
the discrete one as we discussed before. Therefore if one defines a new time t' = Xt, the evolution equation for 
the density matrix acquires an overall factor 1/A. This requires rescaling a in the term involving two Hamiltonians 
by a factor of A to compensate for the fact that the Hamiltonian is rescaled by 1/A. Therefore although the value of 
a has this arbitrariness, the speed at which states decohere is the same, just measured in different time units. 

We still have to ask ourselves how do we connect the time t appearing in the formula for a and the physical time 
one would measure in a laboratory experiment. We notice that there exists a relation between t and iiab of the form 
t = t(yl(tiab)- That is, t is connected with the variable we chose as clock A, and the latter is a physical variable of 
the cosmology that should be connected with the clock in the lab (just like one measures the age of the universe in 
seconds as defined by a regular clock) . This will rescale the left hand side of the evolution equation of the density 
matrix of the form 

(52) 



i9iiab dt 



The prefactor i9t/9tiab will be slowly varying with respect to the time scales of a laboratory experiment. Therefore 
we can take it to be a constant A of the same nature as the one we discussed in the previous paragraph. This leads 
to a definition of the a one would measure in a lab, 

dtuh 

o-iab = —rr—cr- (5o) 
ot 

In the model we are discussing, the closest thing to a "physical" time could be the time of comoving observers. 
This time has an exponential relation to the variable t we have been using. Unfortunately, due to the simplicity of 
this model there is an overall scale ambiguity and the relation between the comoving time and t is defined up to an 
overall constant k. The relation is t = k exp{3^/At / 2) . Therefore the result for the physical a is, 

a > ^.El^ . (54) 

Ak2 exp(3vAt;ah) 

If one were dealing with a universe that is not exactly DeSitter but with a Friedmann model with a power law 
behavior in comoving time, the relation between comoving time and t would not be exponential but be of the form 
^/Xt — (VAtuh)" with a some power, since the Big Bang would have to happen at t = iiab = 0. If that is the case 
and we assume one is making an experiment today, we have that -s/Atiab ^ 1 and therefore (we are setting the speed 
of light c = 1) one has that a ^ tpianck- As expected, the decoherence effect due to the finiteness of the (space)-time 
lattice is of the order of the Planck scale. 
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VI. CONCLUSIONS 

We have shown how the use of a relational time in quantum mechanics leads to a modification of ordinary quantum 
mechanics in which pure states evolve to mixed states. The use of a relational time is now possible in quantum 
gravity due to the use of the consistent discretization framework. One can then proceed to work out estimates of the 
magnitude of fundamental decoherence due to the fact that in a cosmology time has necessarily to be relational. We 
have worked out such estimate in this paper using a simple cosmological model. It is clear that more detailed models 
will be needed to gain confidence in the estimate obtained. 

What are the possibilities of detecting the decoherence due to quantum gravity? We have found that the largest 
amount of decoherence one can expect in a system goes as 

{UJI — a;2)^tpianckilifc (55) 

where the w's are the Bohr frequencies corresponding to the two most separated energy levels of the system and tute 
is how long we wait for the decoherence to appear. For ordinary quantum systems such effect is very small. It is 
interesting to notice however, that a small level of fundamental decoherence has always been desirable in quantum 
mechanics to avoid "Schrodinger cat" type situations, and it had been advocated through ad-hoc proposals like that of 
Ghirardi, Rimini and Weber |1R| | and others. The most promising experiments where one could observe this effect are 
the ones involving macroscopic quantum systems, like Bose-Einstein condensates. We have discussed some of these 
possibilities in reference j2l| , but it appears detectability is not within reach of current technology, largely because the 
systems involve a limited number of atoms. Since the energy spread goes as the number of atoms squared, it might 
be possible that in a few years experiments could reach the desired levels of energy spread in macroscopic quantum 
systems. It remains to be seen if the specific experimental setups will allow the effect to be sufficiently isolated from 
other sources of decoherence due to environmental effects. 

It should be noted that we are claiming that the consistent discretizations allow to solve the problem of time 
in quantum gravity and totally constrained systems in general through the introduction of a relational time. The 
problem of time has many aspects to it. One of these aspects is the identification of a variable that behaves as a 
"good time". In particular, a variable that is "transverse" to the dynamical orbits of the system. Hajfcek et al. [2^ . 
Hartle and Kuchaf Q have emphasized that in many cosmological situations there do not exist time variables 
transverse to the dynamical orbits. Our approach actually addresses this problem. In a nutshell, even in such cases 
the parameter n in the discretization is transverse to the orbits and this is enough to define a correct relational time. 
This in particular implies that our approach solves the "time of arrival" problem in quantum mechanics. We will 
expand on these issues in a forthcoming publication. 

Another aspect of the framework that requires further analysis is the issue of the covariance of the predictions. 
As worked out in this paper, the predictions are particular to a given choice of time. This is unsatisfactory since 
one should expect the theory at least to exhibit Lorentz invariance locally since it is derived from general relativity. 
Establishing this in general relativity is however, delicate, since it requires analyzing situations with spatial (field 
theoretic) degrees of freedom like Gowdy cosmologies which are considerably more involved computationally. We are 
studying these but we do not have results to report at present. An alternative would be to consider model situations, 
for instance, two parameterized relativistic particles with a relativistic interaction, but this also requires further study. 

Summarizing, the use of consistent discretizations of general relativity free the theory from constraints and therefore 
one can introduce a relational time as proposed by Page and Wootters avoiding the main objections to such approach. 
The resulting modified quantum mechanics implies that pure states evolve into mixed states. We have computed 
explicitly the rate for such decoherence. It may play a role in foundations of quantum mechanics, in solving the black 
hole information problem |24j . and may be observed experimentally. 
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